A topological charge selection rule for phase singularities 
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We present an study of the dynamics and decay pattern of phase singularities due to the action 
of a system with a discrete rotational symmetry of finite order. A topological charge conservation 
rule is identified. The role played by the underlying symmetry is emphasized. An effective model 
describing the short range dynamics of the vortex clusters has been designed. A method to engineer 
any desired configuration of clusters of phase singularities is proposed. Its flexibility to create and 
control clusters of vortices is discussed. 
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With the coming of the new century the study and 
understanding of the formation and dynamics of phase 
singularities in nature has become one of the major mile- 
stones in modern physics. The relevance of one particular 
kind of phase singularities, the vortices, has been made 
evident from the variety of scientific fields where they 
play a major role: nuclear physics, condensed- matter 
physics, optics, superfluidity, cosmology... [l| 

This would have not been possible without the ac- 
knowledgement of the important role played by symme- 
try principles in our understanding of nature. With the 
development of physics in the 20th century, symmetry 
considerations evolved from a passive role in which sym- 
metry was a property of the interactions to an active 
role in which symmetry serves to determine the interac- 
tions and the conservation laws themselves, as C.N. Yang 
stated "symmetry dictates interaction " [2j] . 

Although the following results and discussion can be 
applied to any kind of phase singularity we choose one 
particular system to work without any loss of generality, 
the optical vortices. In optics, or in any kind of general 
theory of wave mechanics, vortices are solutions of the 
wave equations characterized by a single phase singular- 
ity, i. e., a screw phase dislocations where the amplitude 
vanishes The phase around the singularity has an 
integer number of windings. For symmetric systems, this 
number, usually denoted as u, is a conserved quantity and 
governs the interactions between vortices as if they were 
endowed with electrostatic charges. Thus, v is referred to 
in the literature as vorticity, topological charge, or wind- 
ing number. These systems have been widely discussed 
in the field of nonlinear optics, both in homogeneous 0] 
and discrete symmetry media [5[. 

Phase singularities in general, and vortex solitons in 
particular, behave and interact with each other in the 
same way as particles do [6|. In fact, skyrmions, soli- 
tons whose topological charge has been identified with 
the baryon number, have been widely used in nuclear 
physics to obtain an effective description of nuclei Q. 
Symmetries and their corresponding conservation laws 



play a crucial role in the description of particles, inter- 
actions and decay patterns, being the Standard Model 
(SM) of electroweak interactions a textbook example 
Therefore, pushing even forward the phase singularity- 
particle analogy, one may wonder whether it is possible 
to characterize the dynamics of the phase singularities as 
particle-like systems by using symmetry arguments. 

In this letter we take advantage of symmetry princi- 
ples and group theory arguments to address the phase 
singularities dynamics and decay pattern, i.e., the for- 
mation of clusters of phase singularities as a result of the 
break-up of the initial one, in the vicinity of an inter- 
face between an homogeneous medium characterized by 
an 0(2) symmetry and a discrete C n one. 

An optical vortex within a medium characterized by 
a continuous 0(2) symmetry presents well-defined an- 
gular momentum i. Such a system being an ideal can- 
didate for new technological developments 0. In sys- 
tems such as optical media with discrete symmetry [Hj or 
Bose-Einstein condensates in 2D periodic traps [10[, the 
0(2) symmetry is replaced by a discrete-point symme- 
try, C n , and the angular momentum becomes ill-defined. 
Nevertheless, although these solutions cannot have well- 
defined angular momentum any longer, certainly all of 
them present neat phase dislocations that can be charac- 
terized by an integer value. In this case one can define an 
integer quantity m, the angular pseudomomentum, which 
is conserved during propagation This behaviour is 
analog to the one of particles crossing into a discrete peri- 
odic medium, usually electrons entering into a crystalline 
system. In these media the linear momentum of the par- 
ticle is no longer a conserved quantity, being replaced on 
this role by the Bloch wavevector [121 ]. 

In the case of homogeneous media, angular momentum 
and vorticity coincide, £ = v. However we shall see that 
this is no longer true whenever the vortex moves into a 
medium characterized by a C n discrete symmetry. To 
discuss its dynamics we will take advantage of discrete 
group theory arguments in the same way it was done in 
Refs. As initial condition we consider the propaga- 
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tion inside an 0(2) medium of a vortex with angular mo- 
mentum t given by <fri(r, (9, z) = e lie fg(r, z), where /g(r, z) 
describes the radial dependence of the vortex amplitude. 
The equation of motion that governs the dynamics of the 
field in both media is the standard nonlinear Schrodinger 
equation given by 



dz 



'2D 



V L (z,y,z)+V NL (\<l>i\)]<l>i (1) 



where V L {x,y,z) = V Cn (r, 0)H(z) + V(r) ( 2 ) [1 " H(z)] 
characterize the media before and after the interface, 
H(x) is the Heaviside step function and Vo(2)( r ) an d 
Vc n (r,0) are rotationally invariant potentials, 0(2) the 
former and C n the latter. V/vl(|0|) is a nonlinear term 
depending on the field modulus. 

In the linear regimen, given an initial condition and 
knowing the eigenstates (eigenmodes), ip p (x,y), and 
eigenvalues (propagation constant), of the system, 
the field fa is completely determined for any value of 
z as fa(r,6,z) = ^2 p e lfJ ' pZ c p (0)ilj p (r,0). In the nonlinear 
case this is no longer true. However, one can perform the 
same expansion in an infinitesimal element dz in which 
the nonlinear term varies slowly and therefore the con- 
clusions holds. The c p (0) coefficients can be obtained 
through the projection of the vortex over the eigenmodes 
of the system. Thus, given the system eigenmodes and 
their projections over the interface the vortex evolution 
is completely determined. 

Group theory allows us to go one step further and ex- 
tract more information about these coefficients. If the 
linear and nonlinear potentials are invariant under the 
C n group, the eigenmodes of the system can be classi- 
fied according to the different group representation. By 
using this, the field over the interface can be rewrit- 
ten as fa(r,6,0) = Y.rn,m C Lm eimeu rnrn(r,6), where m 
numerate the modes with angular pseudomomentum m 
in each representation. Since fa(r, t9, 0) and U m rn(r,0) 
belong to representations of 0(2) and C n respectively, 
they both transform properly under a discrete rotation 
of order n. Thus, by performing the change of vari- 
able — > 6 + 2ir/n one arrives to the symmetry relation 
C e mW = e t2 ^ e - m)/n C e mW . Then, the C e mW coefficients are 
zero unless the condition 



£ = m + kon (fco G Z) 



(2) 



is fulfilled. Even more, in Refs. [ll| it was proved that 
the values of m are further constrained by the order of 
the point-symmetry group C n , being \m\ < n/2 (even n) 
or \m\ < (n — l)/2 (n odd). 

Eq.([2|) also implies that once the angular momentum 
and the symmetry order are fixed (initial conditions of 
our problem), the vortex would only be projected over 
modes belonging to the ra-th representation of the sym- 
metry group, i.e. the dynamics of the vortex will be de- 
scribed in terms of eigenmodes characterized by an an- 
gular pseudomomentum m. This implies that the field 



over the interface can be expressed as a linear combina- 
tion of the form <p mko {r, 9, 0) = J2 w Ct m e m0 U mm (r, 9). 
Umrn{ r , @) are periodic functions in the angular variable 
and therefore they admit a Fourier expansion. Intro- 
ducing this into the definition of (j)mk (r, 0, 0) one arrives 
to <j) mko (r,9,0) = Y. k ,mC k ^e l(m+kn)e Ut m {r). Know- 
ing the expression of the coefficients over the interface 
and summing over the subindex m, the expression of the 
field for any value of z is given by 



</>mk (r,6,z) = ^a^(r, z)e 



i(m-\-kn)0 



(3) 



where a k m (r,z) = ^"'""^CW • Eq.© can 
be understood as a linear combination of vortex func- 
tions whose charges are completely determined by the 
angular pseudomomentum m and the order of the group 
n. The dynamics of each one being determined by the co- 
efficients aj^(r, z). This expression is valid in the linear 
and in the nonlinear case. For the latter, it is impor- 
tant to emphasize that V/vl(|<^|) depends only on the 
modulus of the field, and therefore the eigenmodes have 
well defined symmetry properties under the C n discrete 
group. This makes that Eq.® will be valid within an 
infinitesimal interval dz. If we consider the propagation 
on the interval [0, z] as a sum of subintervals [z^ Zi + dz) 
and we apply the aforementioned arguments to each one 
of them, it can proved that Eq.([3]) would also be valid in 
the nonlinear case. 

The evolution of the radial coefficients is obtained by 
introducing Eq.(j3]) into Eq.(pQ) and integrating out the 
angular part, 



dz 
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kCl n 



(4) 



where V q -k is given by the Fourier expansion of 
Vc n (r,6). The functional form of M q _k will depend on 
the kind of nonlinearity considered. For Kerr nonlineari- 
ties it will take the form M q -k = j(z) J2 P a *k-q+p a p- The 
dynamics described by Eq. ([3]) is completely equivalent 
to the one given by the nonlinear Scrodinger equation. 
However, in this case we have reduced the dimensionality 
of the problem, obtaining a description of the nonlinear 
dynamics in terms only of the radial variables. 

To provide a physical example we consider an opti- 
cal interface separating two 2D dielectric media with 
Kerr nonlinearity, these two media being a homoge- 
neous medium and a 2D square optical lattice. This 
system is equivalent to a 2D BEC in which a pe- 
riodic potential is abruptly switched on. They con- 
stitute an 0(2) — O4 interface described by Eqs. (pQ) 
and © with V NL = |^| 2 , V( r )o(2) = Vb, and 

Vc n (*,v) = ^iE?=oE*p{ ~^ where 
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FIG. 1: Phase calculated for z = 0.30 (a) using Eq. (pQ) and 
(b) Eq. (|4]). The corresponding angular components a^(r, z — 

0. 30) are given in (c). 

(xj, yj) = d (cos (2irj/n, sin(27rj/n)) (d stands for the dis- 
tance between gaussians, and V\ = — 1). 

We consider an initial vortex with charge £ = 5 crossing 
the interface into the C4 symmetric medium. Thus, by 
applying Eq. ([2j) one obtains fco = 1. The initial condition 
in the interface will be given by a^(r, 0) = if k 7^ fco and 
a^(r,0) = fi(r,0) if k = fco, therefore the only nonzero 
amplitude in the interface will be a^. Since we are in- 
terested in the short distances limit (z — > 0,r — > 0), 
Eq. ([3j) can be safely truncated in the vicinity of the 
symmetry axis at second order. Being consistent with 
this approximation we have considered a potential in 
which the Fourier expansion can be truncated in the 
same way, V(r,6) * F_i(r) e -^ + V (r) +V 1 (r)e t40 
where Vo > V±\. We show in Fig. [Tfc) that only a re- 
duced set of angular amplitudes, a^, a^, and a^, play 
a significant role in the description of the dynamics of 
the system in the short distances limit. In Fig. [T] we 
present the phase op of the complex solution </> = \(j)\e llf 
obtained numerically solving Eq. (j4]) and the complete 
non-linear Schrodinger equation Eq. (pQ), to show that 
both simulations provide similar results. From these pic- 
tures it can be observed one singularity lying in the sym- 
metry axis with charge +1, note that according to Eq. (j3]) 
m = = 5) — (fcon = 4) = +1, together with four singu- 
larities with charge +1 symmetrically distributed around 
it. The same breaking pattern is observed for other sim- 
ulations with different £ and n provided that fco = 1, 

1. e., one static singularity in the center with charge m 
and clusters of n singularities with charge +1 moving 
away from the symmetry axis, see EPAPS document 
No. [EPAPS-01.gif] for some examples. Let us estab- 
lish the physical meaning of the angular coefficients a^. 
As shown in Fig. [TJ for r — > and z — > only and 
must be considered. Therefore, in this limit, Eq. ([3]) can 
be written as 0(r, 0, z) « a^(r, z)e ime -f a^(r, ^)e^ m+n ^. 
Note that, if ex r m and oc r m+n this equation 
is transformed in (j)(r,6,z) ~ r m e im6 r m+n e i{m+n9) ^ 
which can be written as <f>(r,9,z) ~ co rn (a + buo n ) with 
co G C, a, 6 G R, and then it correctly describes the 
clusterization pattern observed numerically. The be- 
haviour ex r l m+/cn l is not that surprising since it has 
been demonstrated for stationary solutions, i.e. solutions 
<fi(r,6,z) = e zfiz ilj(r,6) [14]. It has been checked numeri- 




FIG. 2: Phase (a) and dual of the flux field (b) for an initial 
vortex with angular momentum £ — 5 after evolving into a C3 
symmetric medium. 



cally that and a m present this behaviour for r — > 0. 

Once the importance of £ and n has been established, 
let us show the pivotal role played by fco in the mecha- 
nism of the formation of clusters of vortices. To do so 
those solutions where fco > 1 are of particular interest. 
The simplest example will be a £ = 5 vortex entering into 
a C3 symmetric medium, fco = 2. We show in Fig. [2] the 
phase and the dual of the flux field, i.e. that which is 
perpendicular at every point to V(/?, after evolving in a 
Cs medium. From the analysis of the singularities we ob- 
serve one singularity of charge m = — 1 in the center and 
two clusters of vortices, that we name "waves", emanat- 
ing from the center. Each one of these waves is formed by 
a cluster made of three vortices of charge +1. Therefore, 
I fco I stands for the number of waves emanating from the 
center. 

Note that phase dynamics occurs close to the symme- 
try axis, i.e., in the region where the field "feels" the 
presence of the potential. On the contrary, the phase is 
not modified by the existence of the potential in points 
located far away from this region, since the potential 
quickly vanishes at long distances. Consequently, if we 
consider closed curves £ surrounding the symmetry axis 
sufficiently far away from the region where the potential 
is significantly different from zero, the topological charge 
v = 2tt $c ' ^ s ec L ua l t° £. If this curve coincides 
with the boundary of R 2 , Eq. ([2]) can be re-interpreted 
as a total topological charge conservation law. 

This is similar to the case of the total electrical charge, 
whose conservation can be derived from the U(l) local 
gauge invariance of the SM lagrangian. The similarity 
between topological and electrical charges was introduced 
in the analysis of the XY model for interacting spin sys- 
tems in two spatial dimensions [l3j]. Within this model 
the total energy of a set of vortices with total topological 
charge equal to zero is the same as the one obtained for 
a diluted 2D Coulomb gas, therefore identifying a set of 
vortices and antivortices with a neutral plasma formed by 
particles and antiparticles. In an analogous way, Eq. ([2j) 
states that the the decay of an initial high-charged parti- 
cle in a set of charged particles conserves the total charge. 
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FIG. 3: (a) Dual of the flux field for an initial vortex with 
angular momentum £ — 5 after evolving to z — 0.04 into a C4 
symmetric medium tailored as explained in the text, and (b) 
vortices trajectories after evolving to z = 0.1. 



Next we will see that this equation is still satisfied 
if vortex-antivortex pairs are generated. According to 
Eq (|4|) the coupling among the different field angular 
components is given by the potential terms Vk- There- 
fore, by modifying the potential components one can 
manipulate the coupling between the different angular 
components, and hence the clusterization pattern. In 
the aforementioned examples, the potential components 
are decreasing. The linear potential Vl can be engi- 
neered to modify this property, in such a way that the 
clusterization pattern can be controlled. Let us define 
a potential where V2 > V\ with a form Vc n (r, 9) = 
g(r) [Vq + Vi cos(n6>) + V 2 cos(2n<9)]. Within BEC's this 
potential would be obtained by using laser-interference 
based optical traps. In optics, it can be obtained by 
properly controlling the refractive index profile or by in- 
ducing it by means of interfering beams in a photore- 
fractive medium. We show in Fig. [3ja) the dual of the 
flux field of a vortex with angular momentum £ = 5 after 
evolving to z = 0.04 into a C4 symmetric medium char- 
acterized by such a potential. According to the previous 
results, one may expect a static vortex in the center with 
charge +1 and a cluster of four +1 charged vortices em- 
anating away from it. However, the dynamics obtained 
is more involved. One static vortex in the center with 
charge +1 and three waves of vortices have been gen- 
erated. The first two waves are made of four vortices 
with charge +1 each, while the second one is made of 
four —1 charged vortices (antivortices). Therefore, al- 
though the conservation rule holds independently of the 
choice of the potential, the wave dynamics is potential- 
dependent. We plot in Fig. [3](b) the vortices trajecto- 
ries after evolving to z = 0.1 (see EPAPS document No. 
[EPAPS-02.gif] for a 3D animation of the process). It 
can be observed how if the system is allowed to evolve 
beyond a critical distance the four antivortices annihi- 
late with four vortices in the same way particles does, 
radiating their energy away, recovering in this way the 
picture naively expected. This phenomena, well know 
in particle physics, is usually referred to as pair produc- 
tion/annihilation. Quantum number conservation allows 



the creation of pairs of particles-antiparticles providing 
that all the conserved quantum numbers of the produced 
pairs sum zero. Therefore, since according to the Eq. ([2]) 
the total topological charge is conserved, the topological 
charge that goes into the waves (kon) can be distributed 
in any manner consistent with the discrete group sym- 
metry, creating \k v \ + \ko\ waves of n vortices together 
with \k v \ waves of n antivortices. 

The previous results reinforce the particle— phase sin- 
gularity picture. A topological charge conservation rule 
has been identified and analysed. This relation allows to 
describe the vortex decay pattern into a cluster of vortices 
and to fully specify their number and properties, i.e., m 
and the number of vortex, |&o|, and vortex-antivortex, 
I fey |, waves. This process can be engineered by choosing 
the initial angular momentum the symmetry of the 
discrete medium C n , and by tailoring the experimental 
Vc n (r, 6) potential. This allows us to define a control 
mechanism which would allow experimentalist to gener- 
ate and manipulate any desired configuration of waves of 
singularities. Finally, since all the results holds for any 
vortex embedded into a C n symmetric medium, it opens 
the possibility to address a wide variety of phenomena in 
several different physical systems. 
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